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CHANDRA X-RAY SPECTROSCOPY OF THE FOCUSED WIND IN THE CYGNUS X-l SYSTEM 
I. THE NON-DIP SPECTRUM IN THE LOW/HARD STATE 

Manfred Hanke 1,2 , Jorn Wilms 1,2 , Michael A. Nowak 3 , Katja Pottschmidt 4 ' 5,6 , Norbert S. Schulz 3 , Julia C. Lee 7 

ABSTRACT 

We present analyses of a 50 ks observation of the supergiant X-ray binary system CygnusX-l/HDE 226868 
taken with the Chandra High Energy Transmission Grating Spectrometer (HETGS). CygX-1 was in its spec- 
trally hard state and the observation was performed during superior conjunction of the black hole, allowing for 
the spectroscopic analysis of the accreted stellar wind along the line of sight. A significant part of the observa- 
tion covers X-ray dips as commonly observed for Cyg X-l at this orbital phase, however, here we only analyze 
the high count rate non-dip spectrum. The full 0.5-10 keV continuum can be described by a single model con- 
sisting of a disk, a narrow and a relativistically broadened Fe Ka line, and a power law component, which is 
consistent with simultaneous RXTE broad band data. We detect absorption edges from overabundant neutral 
O, Ne and Fe, and absorption line series from highly ionized ions and infer column densities and Doppler 
shifts. With emission lines of He-like Mg XI, we detect two plasma components with velocities and densities 
consistent with the base of the spherical wind and a focused wind. A simple simulation of the photoionization 
zone suggests that large parts of the spherical wind outside of the focused stream are completely ionized, which 
is consistent with the low velocities (<200 km s _1 ) observed in the absorption lines, as the position of absorbers 
in a spherical wind at low projected velocity is well constrained. Our observations provide input for models 
that couple the wind activity of HDE 226868 to the properties of the accretion flow onto the black hole. 

Subject headings: stars: individual (HDE 226868, CygX-1) - stars: winds, outflows - X-rays: binaries - 
accretion, accretion disks - techniques: spectroscopic 


1. INTRODUCTION 

CygnusX-1 was discovered in 1964 (Bowyeretal. 1965) 
and soon identified as a high-mass X-ray binary system 
(HMXB) with an orbital period of 5.6 d (Murdin & Webster 
1971; Webster & Murdin 1972; Bolton 1972). It consists 
of the supergiant 09.7 star HDE 226868 (Walborn 1973; 
Humphreys 1978) and a compact object, which is dy- 
namically constrained to be a black hole (Gies & Bolton 
1982). The detailed spectroscopic analysis of HDE 226868 by 
Herrero et al. (1995) gives a stellar mass M* « 18Mq, lead- 
ing to a mass of Mbh ~10Mq for the black hole, if an in- 
clination i « 35° is assumed. Note that Ziolkowski (2005) 
derives a mass of M* = (40 ± 5 )M 0 from the evolutionary 
state of HDE 226868, corresponding to Mbh = (20 ± 5)Mq, 
while Shaposhnikov & Titarchuk (2007) claim Mbh = (8.7 ± 
0.8 )Mq from X-ray spectral-timing relations. 

CygX-1 is usually found in one of two states which 
are distinguished by the soft X-ray luminosity and spectral 
shape, the timing properties, and the radio flux (see, e.g., 
Pottschmidt et al. 2003; Gleissner et al. 2004a,b; Wilms et al. 
2006): The low/hard state is characterized by a lower lu- 
minosity below lOkeV, a hard Comptonization power law 
spectrum (photon index T ~ 1.7) with a cut-off at high en- 
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ergies (folding energy Lj ol( | ~ 150keV) and strong variabil- 
ity (~30% root mean square, rms). Radio emission is de- 
tected at the ~15mJy level. In the high/soft state, the soft 
X-ray spectrum is dominated by a bright and much less 
variable (only few % rms) thermal disk component, and 
the source is invisible in the radio. Within the classifica- 
tion of Remillard & McClintock (2006), the high/soft state 
of CygX-1 corresponds to the steep power law state rather 
than to the thermal state, as a power law spectrum with pho- 
ton index T ~ 2.5 may extend up to ~ 1 0 MeV (Zhang et al. 
1997; McConnell et al. 2002; Cadolle Bel et al. 2006). Most 
of the time, CygX-1 is found in the hard state, but transi- 
tions to the soft state and back after a few weeks or months 
are common every few years. Transitional or intermediate 
states (Belloni et al. 1996) are often accompanied by radio 
and/or X-ray flares. Similar to a transition to the soft state, the 
spectrum softens during these flares and the variability is re- 
duced. This behavior is called a “failed state transition” if the 
true soft state is not reached (Pottschmidt et al. 2000, 2003). 
Transitional states have occurred more frequently since mid- 
1999 than before (Wilms et al. 2006), which might indicate 
changes in the mass-accretion rate due to a slight expansion 
of HDE 226868 (Karitskaya et al. 2006). 

HMXBs are believed to be powered by accretion from the 
stellar wind. The accretion rate and therefore X-ray lumi- 
nosity and spectral state are thus very sensitive to the wind’s 
detailed properties like velocity, density and ionization. For 
HDE 226868, Gies et al. (2003) found an anti-correlation be- 
tween the Ha equivalent width (an indicator for the wind mass 
loss rate M*) and the X-ray flux. Considering the photoioniza- 
tion of the wind would allow for a self-consistent explanation 
(see, e.g., Blondin 1994): a lower mass loss gives a lower 
wind density and therefore higher degree of ionization due to 
the irradiation of hard X-rays, which reduces the driving force 
of HDE226868’s UV photons on the wind and results in a 
lower wind velocity v, leading finally to a higher accretion 
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rate (oc M*/v 4 , Bondi & Hoyle 1944). However, Gies et al. 
(2008) find suggestions that the photoionization and veloc- 
ity of the wind might be similar during both hard and soft 
state. UV observations allow to probe the photoionization 
in the HDE 226868 /CygX-1 system: Vrtileketal. (2008) re- 
port P Cygni profiles of N V, C IV, and Si IV with weaker ab- 
sorption components at orbital phase <p or h ~ 0.5, i.e., when 
the black hole is in the foreground of the supergiant. This 
reduced absorption, which was already found by Treves et al. 
(1980), is due to the Hatchett & McCray (1977) effect, show- 
ing that those ions become superionized by the X-ray source. 
Gies et al. (2008) model the orbital variations of the UV lines 
assuming that the wind of HDE 226868 is restricted to the 
shadow wind from the shielded side of the stellar surface 
(Blondin 1994); i.e., the Stromgren (1939) zone of CygX-1 
extends to the donor star. This assumption applies, however, 
only to the spherical part of the wind, which might therefore 
hardly contribute to the mass accretion of Cyg X-l. 

As HDE 226868 is close to filling its Roche lobe (Conti 
1978; Gies & Bolton 1986a,b), the wind is not spherically 
symmetric as for isolated stars, but strongly enhanced to- 
wards the black hole (“focused wind”, Friend & Castor 1982). 
The strongest wind absorption lines in the optical are there- 
fore observed at the conjunction phases (Gies et al. 2003). 
Similarly, X-ray absorption dips occur preferentially around 
4>orb = 0, i.e., during superior conjunction of the black hole 
(Balucinska-Church et al. 2000). These dips are probably 
caused by dense, neutral clumps, formed in the focused wind 
where the photoionization is reduced, although recent analy- 
ses also have suggested that part of the dipping activity may 
result from the interaction of the focused wind with the edge 
of the accretion disk (Poutanen et al. 2008). 

The photoionization and dynamics of both the spherical 
and focused wind can also be investigated with the high- 
resolution grating spectrometers of the modern X-ray obser- 
vatories Chandra or XMM-Newton. As none of the previ- 
ously reported observations of Cyg X-l was performed at or- 
bital phase (po T b = 0 and in the hard state, when the wind is 
probably denser and less ionized than in the soft state, the 
Chandra observation presented here allows for the most de- 
tailed investigation of the focused wind to date. 

The remainder of this paper is organized as follows: In §2, 
we describe our observations of CygX-1 with Chandra and 
the Rossi X-ray Timing Explorer ( RXTE ), and how we model 
CCD pile-up for the Chandra-HETGS data. We present our 
investigations in §3: After investigating the light curves, we 
model the non-dip continuum and analyze neutral absorption 
edges and absorption lines from the highly ionized stellar 
wind - and the few emission lines from He-like ions, which 
indicate two plasma components. In §4, we discuss models 
for the stellar wind and the photoionization zone. We summa- 
rize our results after comparing them with those of previous 
Chandra observations of Cyg X-l. 

2. OBSERVATION AND DATA REDUCTION 
2.1. Chandra ACIS-S/HETGS observation 

CygnusX-1 was observed on 2003 April 19 and 20 by the 
Chandra X-ray observatory, see Table 1 . An overview on all 
its instruments is given by the Proposers’ Observatory Guide 
(CXC 2006). In the first four months of 2003, the RXTE All 
Sky Monitor (ASM; see Doty 1994 and Levine et al. 1996) 
showed the source’s 1 .5-12 keV count rate to be generally be- 
low 50 ASM-cps (Fig. 1). At the time of our Chandra obser- 
vation, it was less than 25 cps (0.33 Crab), typically indicative 
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Fig. 1. — The brightness of Cyg X- 1 as seen by the All Sky Monitor (ASM) 
on board RXTE. During the Chandra observation (marked by a line), the 
source was still in its low/hard state: the 1.5—12 keV count rate did not exceed 
25 cps. In spite of the high intrinsic variability, the high/soft state during June, 
July and August can clearly be distinguished, a result which is also found by 
spectral analysis (Wilms et al. 2006). 


TABLE 1 

Observations of Cyg X- 1 


satellite / 

start 

stop 

exposure" 

count rate" 

instrument 

(MJD) 

(MJD) 

(ks) 

(cps) 

Chandra 

52748.70 

52749.28 

(47.2) 

( ) 

MEGil* 

(52748.70) 

(52749.28) 

16.1 

2x27 

HEG±1 C 

(52748.70) 

(52749.28) 

16.1 

2x 17 

RXTE 

52748.08 

52749.18 



PCA d 

52748.74 

52748.78 

3.0 

1456 

HEXTE a+b e 

52748.74 

52748.78 

1.1 

2x186 


“For the Chandra data, the non-dip GTIs (see, e.g., Fig. 4) have been used. 
For RXTE, the 1 1 th orbit was considered. 

6 The Chandra-MEG spectra cover es 0.8-6 keV (1-99 % quantiles). 
r The Chandra-HEG spectra cover Ed 1-7.5 keV (1-99 % quantiles). 
rf The RXTE-PCA data from 4 to 20keV has been used. 
c The RXTE-HEXTE data from 20 to 250 keV has been used. 

of the source being in its low/hard state (Wilms et al. 2006). 

The High Energy Transmission Grating Spectrometer 
(HETGS), containing High and Medium Energy Gratings 
(HEG/MEG; see Canizares et al. 2005) was used to disperse 
X-ray spectra with highest resolution (CXC 2006, Table 8.1): 

AAheg = 5.5 mA and AAMEG=H m A (1) 

As only half of the spectroscopy array of the Advanced CCD 
Imaging Spectrometer (ACIS-S; see Garmire et al. 2003), 
namely a 512 pixel broad subarray, was operated in Timed 
Event (TE) mode, the 6 CCDs could be read out after expo- 
sure times of ff rame = 1 .7 s and the position of each event is 
well determined. Photons from the different gratings can thus 
easily be distinguished due to the different dispersion direc- 
tions of the HEG and the MEG. Even in its low state, how- 
ever, Cyg X-l is so bright that several photons may pile up in 
a CCD-pixel during one readout frame. Both events cannot 
be discriminated and are interpreted as a single photon with 
larger energy. As the undispersed image would have been 
completely piled up, only 10% of those events in a 40 x 38 
pixel window have been transmitted in order to save teleme- 
try capacity. The first order spectra are, however, only mod- 
erately affected, which can be modeled very well (§2.2). The 
alternative to this approach would have been to use the Con- 
tinuous Clocking (CC) mode, where the ACIS chips are read 
out continuously in 2.85 ms, but only the position perpendicu- 
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lar to the readout direction can be determined for every photon 
event. The CC mode was, however, avoided due to difficul- 
ties in the reconstruction of HEG and MEG spectra and other 
calibration issues. 

The undispersed position of the source is required for the 
wavelength calibration of the spectra. We redetermined it 
to RA = 19 h 58 m 2 1 ?67, 6 = +35° 12'5."83 from the intersec- 
tion of HEG and MEG arm and the readout streak (Ishibashi 
2006). Afterwards, the event lists were reduced using the 
standard software from the Chandra X-ray Center (CXC), 
CIA0 3.3 with CALDB 3. 2. 3. 8 Exceptionally narrow extrac- 
tion regions had to be chosen as the background spectrum 
would otherwise have been dominated by the dispersed ex- 
tended X-ray scattering halo around the source (Xiang et al. 
2005). The further analysis was performed with the Interac- 
tive Spectral Interpretation System ISIS 1 .4.9 (Houck 2002). 9 

We use the four first order MEG and HEG spectra (with two 
dispersion directions each, called +1 and -1 in the following) 
which provide the best signal-to-noise ( S/N ) ratio. The “sec- 
ond order spectra” are dominantly formed by piled first order 
events which reach the other order sorting window of data 
extraction (defined in dispersion-energy space) when the en- 
ergy of two first order photons accumulates. This effect is 
most evident for the MEG, whose even dispersion orders are 
suppressed by construction of the grating bars (CXC 2006). 

2.2. Model for pile-up in grating observations 

For the first order spectra, pile-up causes a pure reduction of 
count rate: A multiple event, i.e., the detection of more than 
one photon in a CCD pixel during one readout time, which 
cannot be separated, is either rejected by grade selection dur- 
ing the data-processing or migrates to a higher order spec- 
trum. The Poisson probability for single events in a 3 x 3 pixel 
event-detection cell i (see Davis 2002, 2003; CXC 2005) is 

Pfi) = A(f) • exp(-A(/)) , (2) 

where the expected number of events, A (i) = 70 • Cj„,(;j, is 
given by the total spectral count rate, C to t(f), at this position 
(in units of counts per A and s), and where the constant 70 is 

70 = 3AA-f frame (3) 

where A A is the resolution of the spectrometer of Eq. (1), and 
fframe is the frame time. We therefore describe the pile-up in 
the first order spectra with the nonlinear convolution model 
simple_gpile2 in ISIS, 10 which exponentially reduces 
the predicted count rate C(A) according to Eq. (2): 

C'(A) = C(A) •exp(- 7-C tot (A)) (4) 

Here, the scale 7 « 70 is left as fit-parameter, and C to t(A) takes 
also the photons into account which are dispersed in a higher 
order m < 3. The count rates are estimated from the corre- 
sponding effective areas A,„ and the assumed photon flux S: 

Ctot(A) = V 3 AJX/ni) ■ S(X/m) (5) 

z — J m= 1 

simple_gpile2 is based on the simple_gpile model 
(CXC 2005; Nowak et al. 2008), which parameterizes the 
strength of pile-up by the (maximum) pile-up fraction 
p= 1 - exp (—7 • max{Ctot}), while using the parameter 7 of 

8 See http://cxc.harvard.edu/ciao3. 3/. 

9 See http://space.mit.edu/cxc/isis/. 

10 The ISIS/S-LANG code for simple_gp±le2 is available online at 
http://pulsar.sternwarte.uni-erlangen.de/hanke/X-ray/code/simple_gpile2.sl . 



Fig. 2. — Pile-up in the HEG (black curves) and MEG (gray curves) data: 
The lower spectra, which are shifted by a factor of 0.4, show the uncorrected 
data. The dashed line shows the model (free of pile-up). The MEG spectrum 
suffers from significant pile-up losses around 2 keV, where the highest count 
rate is obtained. The upper spectra show the pile-up corrected spectra. Note 
that we show ISIS’ (model independent) flux spectra only for illustration; 
simple_gpile(2) operates on the count rates predicted by a model. 

simple_gpile2 avoids to have a non-loacl model which 
depends on the flux at the position of highest pile-up. 

The effect of pile-up is stronger in the MEG spectra than 
in the HEG spectra due to the lower dispersion and higher 
effective area of the MEG. The apparent flux reduction is 
most significant near 2keV where the spectrometer has the 
largest efficiency and the highest count rates are obtained 
(see Fig. 2). It can, however, excellently be modeled with 
simple_gpile2. When fitting a spectrum, there is always 
a strong correlation between the pile-up scale 7 and the cor- 
responding flux normalization factor, e.g., the relative cross 
calibration factor c introduced in §3.2. The best-fit values for 
the pile-up scales 7 found in our data analysis (see Table 2) 
are only slightly larger than expected from Eq. (3) (namely 
7o,meg±i = 5.6 x 1 CT 2 s A and 7o,heg±i = 2.8 x 10“ 2 s A), and 
the calibration factors c are consistent with 1, except for the 
HEG-1 spectrum, for which both the largest 7—70 and c were 
found. Given the presence of the 7 -c- correlation, we consider 
the latter to be a numerical artefact. 

According to the simple_gpile2 model, the MEG+1 
spectra suffer from >30 % pile-up for 6 A < A < 8.3 A, peak- 
ing at pmeg+i = 45% in the Si XIII f emission line at 6.74 A. 
Except for some emission lines, among them the Fe Ka line, 
the continuum pile-up fraction of the HEG spectra is below 
17 %. For the HEG+1 spectrum, the reduction is <10% out- 
side the ranges 2.09 A < A < 4.08 A and 6.05 A < A < 6.93 A. 


2.3. RXTE observation 

While Chandra's HETGS provides a high spectral resolu- 
tion, the energy range covered is rather limited. Within the 
framework of our RXTE monitoring campaign, the broad 
band spectrum of CygX-1 was measured regularly, i.e., at 
least biweekly, since 1998 (see Wilms et al. 2006, and refer- 
ences therein). The observation on 2004 April 19 was ex- 
tended to provide hard X-ray data simultaneously with the 
Chandra observation. More than one day was covered by 17 
RXTE orbits of ~47 min good time, interrupted by ~49 min 
intervals when CygX-1 was not observable due to Earth oc- 
cultations or passages through the South Atlantic Anomaly 
(SAA), see Fig. 3. 
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Fig. 3. — Coverage of the simultaneous RXTE and Chandra observations. The plot shows the background- subtracted light curves with a time resolution of 64 s 
on a separate y-axis each. Top: RXTE - HEXTE (20-250 keV). The count rates have been corrected for the detector dead time. Center: RXTE -PC A (4—20 keV), 
normalized by the number of active PCUs. Bottom: Chandra - HETGS (0.5-12 keV, only first order events; see Fig. 4 for more details). The numbers label the 
RXTE orbits, vertical lines mark the first part of the non-dip spectrum (see §3.1 and Fig. 4), which completely covers RXTE orbit 11. 
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Fig. 4. — Top: Full 0.5-12 keV band Chandra light curve. Bottom: Ratio 
of 0.7-1 keV band and 2.1-7.2keV band count rates. Absoiption dips - at 
first compact, then with complex substructure - show up with a reduced flux 
and spectral hardening. Count rates >82.7 cps define the non-dip data (dark). 


Fig. 5. — The RXTE 4-250 keV broad band continuum spectrum, as mea- 
sured with PC A below 20keV and above 20keV with HEXTE, can be de- 
scribed by a broken power law with high energy cutoff and a weak iron line, 
see Table 2. The HEXTE spectra shown are renormalized to match the PC A 
flux, as the absolute calibrations of those two instruments differ by ~15%. 


Data in the 4-20 keV range from the Proportional Counter 
Array (PCA; Jahoda et al. 1996) and in the 20-250 keV range 
from the High Energy X-ray Timing Experiment (HEXTE; 
Gruber et al. 1996) were used. The data were extracted using 
HEASOFT 6. 3.1, 11 following standard data screening proce- 
dures as recommended by the RXTE Guest Observer Facility. 
Data were only used if taken more than 30 minutes away from 
the SAA. For the PCA, only data taken in the top-layer of 
the proportional counter were included in the final spectrum, 
and no additional systematic error was added to the spectrum. 
During the observation, different sets of Proportional Counter 
Units (PCUs) were operative. During the 11 th orbit, exten- 
sively used in this work, PCUs 1 and 4 were off. 

3. ANALYSIS 
3.1. Light cur\’e 

The Chandra observation covers a phase range between 
~0.93 and ~0.03 in the 5.599829 d binary orbit (Gies et al. 
2003, whose epoch is HJD 2451730.449 ± 0.008). The top 
panel of Fig. 4 shows the light curve of first order events 
(MEG±1, HEG±1) in the full band accessible with Chan- 
dra -HETGS. During several dip-events, the flux is consider- 
ably reduced. The absorption dips distinguish themselves also 
by spectral hardening (bottom panel of Fig. 4). We extract 
a 16.1 ks non-dip spectrum, which is the subject of this pa- 
per, from all times when the total count rate exceeds 82.7 cps; 
these are indicated by dark points in Fig. 4. The analysis of 
dip-spectra will be described in a subsequent paper. 

11 See http://heasarc. gsfc.nasa.gov/lheasoft/. 


The RXTE light curve shows considerable variability as 
well. Although the dips are more obviously detected with 
Chandra in the soft X-ray band, similar structures are also 
seen with RXTE -PCA or even -HEXTE, especially in the last 
RXTE orbits of these observations (Fig. 3). We chose to infer 
the non-dip broad band spectrum from the RXTE data taken 
during the 1 1 th orbit, which was performed entirely during the 
first part of the non-dip phase. Other parts are interrupted by 
dips, occultations, or have nonuniform PCU configuration. 

3.2. Continuum spectrum 

The broad band spectrum of Cyg X-l in the hard state can 
be described by a broken power law with exponential cut- 
off. Since the parameters of this phenomenological model 
are correlated with those of physical Comptonization models 
(see, e.g., Wilms et al. 2006, Fig. 1 1), we are justified in using 
the aforementioned simple continuum model for this paper fo- 
cussing on spectroscopy of the wind. We describe the whole 
0.5-250 keV spectrum consistently with one broken power 
law model, i.e., there is no need to fit the continuum locally. 

Taking pile-up into account (§2.2), the non-dip Chandra 
spectra alone can already be described quite well by a weakly 
absorbed, relatively flat power law spectrum with a pho- 
ton index Ti = 1.51 ±0.01 (Table 2). This result is consis- 
tent with the fact that the break energy of the broad band 
broken power law spectrum is found at Ebreak = 9keV, i.e., 
the Chandra data are virtually entirely in the regime of the 
(steeper) photon index T i . The onset of the exponential cut- 
off (with folding energy Efoid) is at E cut = 24keV and thus 
also well above the spectral range of Chandra. As it is 
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TABLE 2 

Fit parameters of the continuum non-dip hard state spectrum 
ofCygX-1 




fit to the 

joint fit to both the 

parameter 

unit 

Chandra 

Chandra and RXTE 



spectra only 

spectra 

photoabsorption 




Ah 

10 21 cm' 2 

3.52 ±0.04 

5.4 ±0.4“ 

(broken) power law 



G (HETGS) 


1.51 ±0.01 

1.60±0.01 

r, (PCA) 



1.73 ± 0.01 

^break 

keV 


9.°T?;| 

r 2 



1.50±0.01 

norm 

s -1 cm -2 keV -1 

1.15 ±0.01 

1.33±0.03 

high energy cut-off 



£cut 

keV 


24f 2 

^fold 

keV 


204 ±9 

disk black body 
4 disk (Eq. 6) 

10 3 


23+ 17 

-12 

kTcoi 

keV 


0 ?5 +0 03 
u - ZJ -0.02 

narrow iron Ka line 



^0, narrow 

keV 


6 400 1+ 0 0004 
o.h-uui_q 00 g2 

^narrow 

eV (!) 


<0.1 

-^narrow 

1(T 3 s -1 cm' 2 


1.0±0.3 

broad iron Ka line 



^0, broad 

keV 


6.4 (fixed) 

Abroad 

keV 


0.6 ±0.1 

Abroad 

10' 3 s' 1 cm' 2 


4 3+2.7 

relative flux calibration (constant factor) 


CMEG-I 


1 (fixed) 

1 (fixed) 

cmeg+i 


1.00 ±0.01 

0.99 ±0.01 

cheg-i 


1.04 ±0.01 

1.04±0.01 

C HEG+1 


1.01 ±0.01 

1.02±0.01 

CPCA 



1.18 ± 0.03 

CHEXTE 



1 03 +0 ' 02 
1 - u - 0.03 

pile-up scales 




7MEG-1 

10“ 2 s A 

5.6 ±0.1 

5.6±0.1 

7MEG+I 

10' 2 s A 

6.0 ±0.1 

5.9±0.1 

7HEG-1 

10' 2 s A 

4.5 ±0.3 

4.6 ±0.3 

7HEG+1 

10' 2 s A 

3.5 ±0.3 

3.8±0.3 

fit-statistics 




x 2 


11745 

12180 

dof 


1 1274 fe — 10 c 

11274 i +293 6 -25 c 



1.04 

1.06 

absorbed flux (pile-up corrected) 



So.5-10keV 

photons s 1 cm 2 


1.4 

To.5-10keV 

10' 9 erg s' 1 cm' 2 


7.4 

So.5-250keV 

photons s -1 cm -2 


1.8 

To. 5-250 keV 

10 -9 erg s -1 cm -2 


24 

unabsorbed luminosity, assuming d 

= 2.5 kpc (Ninkov et al. 1987) 

Eo.5-10keV 

10 37 erg s 1 

0.67 

0.78 

Eq. 5-250 keV 

10 37 erg s' 1 


3.1 

Error bars indicate 90% confidence intervals for one interesting parameter. 

"In the joint fit, photoabsorption was described with tbnew model 
(Wilms et al. 2000; Juett et al. 2006b) and the best-fitting abundances (see 
§3.3). 

b All data have been rebinned to contain > 50 counts bin -1 . 

c The model contains actually more parameters, as 

the absorption lines have 


already been included (see Tables 4—6). 



Fig. 6. — The C/iandra-HEG spectrum in the Fe line region. The model in- 
cludes both the narrow and the broad Ka emission line, the latter as required 
by the simultaneous RXTE - PCA data, and the Fe XXV r absorption line. 

known that there are cross-calibration uncertainties between 
Chandra and RXTE (Kirsch et al. 2005), we use constant 
factors c, for the relative flux calibration of every spectrum 
and also separate parameters T i (HETGS) and T i (PCA), for 
which we find similar values within the joint model (see Ta- 
ble 2). In order to describe a weak soft excess, we add a 
thermal disk component, which accounts for --9 % of the un- 
absorbed 0.5-10 keV flux. The disk has a color temperature 
of kT co i = 0.25 keV; similar to that found by Makishima et al. 
(2008) and Baluciriska-Church et al. (1995), who described 
the soft excess with a kT = 0.13 ±0.02keV blackbody only, 
but the temperature 7' co | may be too high by a factor f > 1.7 
(Shimura & Takahara 1995). The norm parameter of the 
diskbb model is 


Adisk = diskbb . norm 


/ Z^coi/km \ 

\t//10kpc ) 


■ cos 9 , 


(6) 


where d is the distance and 6 is the inclination of the 
disk, which can deviate from the orbital inclination i ss 35° 
(Herrero et al. 1995) as the disk may be precessing with a tilt 
S ss 37° (Brocksopp et al. 1999). R co i is related to the inner 
radius of the disk, R m = r]g(i)f 2 R co \ (with rj ss 0.6-0. 7 and 
g(i) ss 0.7 -0.8; Merloni et ai. 2000). In spite of these un- 
certainties and the large statistical error of ,4 disk (more than 
50 %) the inner disk radius can be estimated to 1 -5Rs < Rm < 
10/?s if a distance d « 2.5kpc (Ninkov et al. 1987) and a 
Schwarzschild radius R$ « 30km (Herrero et al. 1995) are as- 
sumed. Thus, the disk is consistent with extending close to the 
innermost stable circular orbit (ISCO). 

The good S/N ratio afforded by the RXTE - PCA data 
clearly reveals an iron fluorescence K a line. While the instru- 
mental response of the proportional counters does not allow 
for a resolution of the line profile details, i.e., whether it is 
narrow or relativistically broadened, the Chandra-HEG spec- 
tra (Fig. 6) do resolve a strong narrow component at 6.4 keV. 
Nevertheless, since the integrated flux of the Chandra mea- 
sured line is insufficient to account for all of the PCA residu- 
als in this region, we include an additional broad feature to our 
modeling, which is also compatible with the Chandra spec- 
trum. Given the relatively low S/N ratio at these energies, we 
do not model the broad iron line with a proper physical model 
such as a relativistically broadened line, but use a Gaussian 
with its energy fixed at 6.4 keV, as the latter is hardly con- 
strained by the data. These results illustrate the synergy of the 
simultaneous observation with complementary instruments, 
as the combination of narrow and broad line could only be 
revealed by the analysis with a joint model. 
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Wavelength [A] 

Fig. 7. — The spectral region around the Fe L-edges. The top panel shows 
the model flux in units of 10 -5 ph. s -1 cm -2 A -1 (solid line). The dashed lines 
describe the unshifted model, see text. The second panel shows the count 
rate of all 4 HETGS spectra (MEGdil, HEGzbl), which have been combined 
with a resolution of 10 mA, and the folded model. The bottom panel displays 
the residuals x = (data-model)/error. 



Fig. 8. — As in Fig. 7, but for the spectral region around the O K-edge. 
Both MEGd= 1 spectra have been combined with a resolution of 0. 1 A. 

Our global model for the continuum spectrum now enables 
us to address the features of the high-resolution Chandra spec- 
tra, which is the topic in the remainder of this section. 

3 . 3 . Neutral absorption 

Absorbing columns can be measured most accurately from 
the discrete edges in high-resolution spectra at the ionization 
thresholds. We detect the most prominent L-shell absorp- 
tion edge of iron and the K-shell absorption edges of oxy- 
gen and neon. Juett et al. (2004, 2006a) have inferred the fine 
structure at those edges from earlier Chandra-HETGS obser- 
vations of CygX-1 and other bright X-ray binaries: The Fe 
L 2 and F 3 edges, due to the ionization of a 2p , / 2 or a 2p 3 / 2 
electron, respectively, are separately detectable at 17.2 A and 
17.5 A. The O K-edge at 22.8 A is accompanied by (ls-^2p) 
Ka and higher (Is— > up) resonance absorption lines. The K a 
line occurs at 23.5 A for neutral O I and at lower wavelengths 
for ionized oxygen. In the case of neon, neutral atoms have 
closed L-shells, such that Nel only shows a (Is— >3p) K/3 
absorption line close to the K-edge at 14.3 A, while ionized 
neon also shows Ka absorption lines, e.g., Ne II at 14.6 A and 
Neill at 14.5 A. Improved modeling of the neutral absorp- 
tion that takes these features into account has recently been 



14.2 14.4 14.6 

Wavelength [A] 


Fig. 9. — As in Fig. 7, but for the spectral region around the Ne K-edge. 
All spectra have been combined with a resolution of 5 mA. The dotted lines 
show the absorbed continuum model without additional absorption lines. The 
light dashed line shows the model without the absorption of Ne. 

included in the photoabsorption model tbnew 12 (Juett et al. 
2006b), an extension of the commonly used tbvarabs 
model (Wilms et al. 2000). 

As part of the spectral model for the whole continuum, the 
tbnew model can be used to describe the absorption edges 
detected with the Chandra observation of CygX-1 discussed 
in this paper. Figure 7 shows the Fe F-edges requiring a 
blueshift by AA/Ao • c = (540 ± 230)kms _1 of the tbnew 
model, which relies on the cross-sections of metallic iron 
measured by Kortright & Kim (2000). Unlike Schulz et al. 
(2002), Miller et al. (2005) and Juett et al. (2006a) have also 
found that the Fe F 3 edge requires a small shift; their mean po- 
sition of maximum optical depth is Afcl, = 17.498 A, but our 
value Ap e L 3 = (17.469 ±0.014) A is still lower. The shift could 
be caused by the Doppler effect due to a moving absorber, 
by a modified ionization threshold due to chemical bonds, or 
ionization of the iron atoms (van Aken & Fiebscher 2002). In 
an analysis of the CygX-1 high/soft and low/hard state, fo- 
cused on this spectral region, Fee et al. (2008) find that the 
Fe F-edges here can likely be modeled by a heterogeneous 
combination of gas and condensed matter of iron in combina- 
tion with oxygen local to the source environment. If, as sug- 
gested by these authors, the magnitude of the shift is due to 
molecules and/or dust, this shift is one identifying signature 
of the composition and charge state of the condensed state 
material. Such direct Chandra X-ray spectroscopic detection 
of dust via its associated edge structure was first suggested 
for observations of the Fe F-edge in the AGN MCG-6-30-15 
(Fee et al. 2001) and for the observations of the Si K-edge in 
the microquasar GRS 1915+105 (Fee et al. 2002). The latter 
study associated the observed Si K-edge structure with Si0 2 , 
although the origin - source environment or Chandra CCD 
gate structure - was unclear in this case. 

The blueshift of (100 ± 130) kms^ 1 which has been deter- 
mined for all other absorption edges is consistent with zero. 
The O K-edge can only be seen in the MEG spectra after heav- 
ily rebinning (Fig. 8). Nevertheless, the Ka resonance ab- 
sorption line of O I is clearly detected. The region around the 
Ne edge (Fig. 9) is dominated by Fe XVIII absorption lines, 
possibly blending with the Ka absorption lines of Ne II and 
Ne III. No other strong edges are clearly visible in the spec- 

12 The code for tbnew is available online at http://pulsar.sternwarte.uni- 
erlangen.de/wilms/research/tbabs/ for beta testing and will soon be released. 
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TABLE 3 

Column density N and abundance A = N/Nh of neutral absorbers detected along the line of sight towards CygX-1 







— analysis by — 





this work 


Schulz et al. (2002) 

Juett et al. (2004, 2006a) 




observation (obs.) # 3814 

obs. # 107 obs. # 107 

obs. # 3407 

obs. # 3724 



(G).5-10kev = 0.7 X 10 37 

ergs *) 

(io.5-iOkev = 1.6x 10 37 ergs' 1 ) 

(3.1 X 10 37 erg s' 1 ) 

(soft state) 

element 

12 + logA>SM e „ t » 

^element /^Fe 

A , /aISM b 

^element /^element 

^element C 

^element ^element 

^element 

^element ^ 





( 10 17 crrT 2 ) 

(10 17 cm' 2 ) (10 17 cm' 2 ) 

(10 17 cm~ 2 ) 

(10 17 cm' 2 ) 

0 

8.69 

15.4±0.3 

1.45±0.01 

38-1T84 

39.2 ±2.3 63 ±14 

24 ±3 

29 ±3 

Ne 

7.94 

3.91 ±0.18 

2 07+ 004 
z ' u/ -0.03 

9.6 ±0.2 

9.43 ±0.32 7.1 ±1.0 

7 4 + 0.7 
' ' -0.3 

8.6 ±0.7 

Na 

6.16 

1.26 ±0.15 

40 ± 2 

3.1 ±0.2 




Mg 

7.40 

2 . 1 6 ± 0 . 1 3 

4.0 ±0.1 

5.3±0.2 

3.7 ± 1.3 



Al 

6.33 

<0.13 

< 1.8 

<0.2 




Si 

7.27 

1.45 ±0.17 

3.8 ±0.3 

3.8±0.3 

2.3 ± 1.8 



S 

7.09 

0.6 ±0.2 

2 4+°- 5 
-0.3 

1 6 +0 - 3 
1 • - 0.2 




Ar 

6.41 

<0.06 

< 1.0 

<0.1 




Ca 

6.20 

<0.05 

< 1.3 

<0.1 




Cr 

5.51 

<0.18 

1 3 +2 - 7 
1 -1.3 

0 02 +0 05 
- 0.02 




Fe 

7.43 

1 

1.75 ±0.03 

2.52 ±0.04 

1.5 ±0.1 l-47jj>; 23 

1 Ob 4 ^ 07 
1 - UD -0.11 

0.96 ±0.09 


“We use elemental abundances in the interstellar medium, A*™ ent according to Wilms et al. (2000), i.e., xspec_abund ( "wilm" ) ; in ISIS. 

6 The relative elemental abundance, A e i ement /A*™ c t , is a parameter of the tbnew model. 

^element is the product of Aeiemem and Nh (see Tab?e 2). 

^Schulz et al. (2002) use the cross sections from Vemer et al. (1993), and Kortright & Kim (2000) for Fe. 

“duett et al. (2004, 2006a) use the cross section of Gorczyca & McLaughlin (2000) for O, Gorczyca (2000) for Ne, and Kortright & Kim (2000) for Fe. 


trum. The Na K-edge (at 11.5 A; Verner& Yakovlev 1995) 
blends with an absorption line due to the 2s — > 3 p excitation 
of Fe XXII. The Mg K-edge (at 9.5 A) blends with the Ne X 
Ly S absorption line. The Si K-edge (at 6.7 A) is strongly af- 
fected by pile-up and blends with the Mg XII Ly f3 absorption 
line. The S K-edge (at 5.0 A) is relatively weak. Neutral ab- 
sorption from these elements is nevertheless required within 
the tbnew model. 

The results for the individual abundances. A, = A; /Ah, and 
resulting column densities, A, are presented in Table 3. The 
alpha process elements O, Ne, Mg, Si, S are overabundant 
with respect to the interstellar medium (ISM) abundances as 
summarized by Wilms et al. (2000) and therefore suggest an 
origin in the system itself. The total column densities are also 
compared with the values obtained by Schulz et al. (2002) and 
Juett et al. (2004, 2006a) from other Chandra observations of 
CygX-1. Table 3 includes the corresponding source lumi- 
nosities if they are reported in the literature (see also §4.4 and 
Fig. 15). The X-ray flux was highest during the soft state 
observation # 3724. The column densities confirm the conjec- 
ture of Juett et al. (2004) that a higher (soft) X-ray flux ionizes 
material local to the CygX-1 system and reduces the neutral 
abundances. 

The inferred hydrogen column density Ah (see Table 2) 
is in very good agreement with that from ASCA obser- 
vations during the soft state in 1996, namely Ah = (5.3 ± 
0.2) x 10 21 cm -2 (Dotani et al. 1997) and also with Ah = 6.2 x 
10 21 cm -2 obtained from two other different Chandra obser- 
vations (Schulz et al. 2002; Miller et al. 2002, see also §4.4). 
We note that the large column density towards Cyg X-l found 
by many online tools, (7. 2-7. 8) x 10 21 cm -2 , is obtained from 
a coarse grid - with (0.675 deg) 2 pixel size - of Ah mea- 
surements at 21 cm (Kalberla et al. 2005), which does not re- 
solve the strong variations of A H in the region around Cyg X- 1 
(Russell et al. 2007). 


3.4. Absorption lines of H- and He-like ions 
The high-resolution spectra reveal a large number of ab- 
sorption lines of highly ionized ions. The 1.5- 15 A range 
is shown in Fig. 10 as the ratio between the data and the 
continuum-model. As the line profiles are not fully resolved, 
we model each line / with a Gaussian profile G/. In terms of 
the continuum flux model, A cont , the global model reads 


F( A) = e- T(A > ■ F cont (A) 


1±^G,(A) 

i 


’ Fcont(A) 


(7) 


From a Gaussian’s centroid wavelength A and the rest wave- 
length Ao of the identified line, the radial velocity 

v = (A-Ao)/A 0 -c (8) 


of the corresponding absorber can be inferred. With the defi- 
nition in Eq. (7), the norm of Gi is just the equivalent width: 


W XJ := J[l-F,(X)/F C 0 Dt (X)] dX = j G,( A)dA (9) 


W\.i is related to the absorber’s column density, as a bound- 
bound transition i — > j (with rest frequency ;z 0 and oscillator 
strength ff) in an absorbing plasma with column density A, 
creates the following line profile (see Mihalas 1978, §9-2): 


e ~rM 


exp < -A, 


fjj 

m e c Aizq V 47t Azz d ’ 


H 


v-v Q \ 1 

A u D )j 


(10) 


Assuming pure radiation damping, the damping constant F 
equals the Einstein coefficient A j,. The Doppler-broadening 
Au d = v 0 • Co /c is given by Co 2 = 2 kT / m Itm +v^ urb , i.e., is due to 
the thermal and the turbulent velocity of the plasma. For opti- 
cally thin lines (with r( v) <C 1), the equivalent width is inde- 
pendent of Avd, such that the absorbing column density can 
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TABLE 4 

Overview on the detected lines from H- and He-like ions: theoretical rest wavelengths in A 



transition 

hydrogen-like (1 electron) 

O 

VIII 

Ne 

X 

Na 

XI 

Mg 

XII 

Al 

XIII 

Si 

XIV 

s 

XVI 

Ar 

XVIII 

Ca 

XX 

Fe 

XXVI 

Ni 

XXVIII 

Ly a 

1s( 2 S 1/2 ) 

-t- 

2 P(” P 3/2,l/2) 

18.97 

12.13 

10.03 

8.42 

7.17 

6.18 

4.73 

3.73 

3.02 

(1.78) 

1.53 

Ly /3 

1s( 2 S 1/2 ) 

-> 

3p ( 2 P3/2,l/2) 

16.01 

10.24 

8.46 

7.11 

6.05 

5.22 

3.99 

3.15 

(2.55) 

1.50 

1.29 

Ly 7 

1s( 2 S 1/2 ) 

-> 

4 P( 2p 3/2,l/2) 

15.18 

9.71 

8.02 

(6.74) 

5.74 

4.95 

3.78 

2.99 

2.42 

1.43 

1.23 

Ly S 

Is ( 2 S 1/2 ) 

-> 

5p( 2p 3/2,l/2) 

14.82 

9.48 

7.83 

6.58 

(5.60) 

(4.83) 

3.70 

(2.92) 

(2.36) 

(1.39) 

( 1 . 20 ) 

Lye 

1s( 2 S 1/2 ) 

-> 

6p( 2p 3/2,l/2) 

(14.63) 

9.36 

( 7.73) 

6.50 

(5.53) 

(4.77) 

(3.65) 

( 2 . 88 ) 

(2.33) 

(1.37) 


LyC 

Is ( 2 S 1/2 ) 

-> 

7 P ( 2p 3/2,l/2) 

(14.52) 

9.29 

( 7.68) 

6.45 

(5.49) 

(4.73) 

(3.62) 

( 2 . 86 ) 

(2.31) 

(1.36) 


Lyr? 

1s( 2 S 1/2 ) 


8 P ( 2p 3/2,l/2) 

(14.45) 

9.25 

( 7.64) 

(6.42) 

(5.47) 

(4.71) 

(3.60) 

(2.85) 

(2.30) 

(1.36) 


Ly 9 

1s( 2 S 1/2 ) 

-+ 

9 P ( 2p 3/2,l/2) 

(14.41) 

( 9.22) 

(7.61) 

(6.40) 

(5.45) 

(4.70) 

(3.59) 

(2.84) 

(2.29) 

(1.35) 


limit 

1s( 2 S 1/2 ) 

-► 

oo 

(14.23) 

(9.10) 

( 7.52) 

6.32 

5.38 

(4.64) 

(3.55) 

(2.80) 

(2.27) 

1.34 

(1.15) 



transition 

0 

Ne 

Na 

Mg 

Al 

Si 

S 

Ar 

Ca 

Fe 

Ni 


helium-like (2 electrons) 

VII 

IX 

X 

XI 

XII 

XIII 

XV 

XVII 

XIX 

XXV 

XXVII 

f [em.] 

Is 2 (‘So)*- 

- 1s2s( 3 Si) 

22.10 

(13.70) 

11.19 

9.31 

7.87 

6.74 

5.10 

(3.99) 

(3.21) 

(1.87) 


i [em.] 

Is 2 (‘So)*- 

- ls2p( 3 P u ) 

21.80 

(13.55) 

11.08 

9.23 

7.81 

(6.69) 

5.07 

3.97 

3.19 

( 1 . 86 ) 

1.60 

r = He a 

Is 2 (‘S 0 )- 

* ls 2 p( lp 1 ) 

21.60 

13.45 

11.00 

9.17 

7.76 

6.65 

5.04 

3.95 

3.18 

1.85 

1.59 

He/3 

Is 2 (‘S 0 )- 

* ls3p(‘ p i) 

18.63 

11.54 

9.43 

7.85 

6.64 

5.68 

4.30 

3.37 

2.71 

1.57 

(1-35) 

He 7 

Is 2 (‘So)- 

* ls4p(‘ p i) 

(17.77) 

11.00 

8.98 

7.47 

6.31 

5.40 

4.09 

3.20 

(2.57) 

1.50 

(1.28) 

He S 

Is 2 (‘So)- 

* ls5p(‘ p i) 

(17.40) 

10.77 

8.79 

7.31 

(6.18) 

(5.29) 

4.00 

(3.13) 

(2.51) 

1.46 

(1.25) 

Hee 

Is 2 (‘So)- 

-* ls 6 p(‘Pi) 

(17.20) 

10.64 

( 8.69) 

(7.22) 

( 6 . 10 ) 

5.22 

3.95 

(3.10) 




HeC 

Is 2 (‘So)- 

* ls7p(‘P]) 

(17.09) 

10.56 

( 8.63) 

(7.17) 

(6.06) 

(5.19) 

(3.92) 





Her; 

Is 2 (‘So)- 

* ls 8 p(‘P 1 ) 

(17.01) 

(10.51) 

( 8.59) 

(7.14) 

(6.03) 

(5.16) 

(3.90) 





limit 

Is 2 (‘S 0 )- 

+ Is oo 

16.77 

10.37 

8.46 

(7.04) 

5.94 

5.09 

(3.85) 

(3.01) 

2.42 

1.40 

( 1 . 20 ) 


Lines with (wavelengths in parentheses) are not detected in our Chandra-HETGS observation of CygX-1, while lines indicated with bold wavelengths are 
clearly detected and those with underlined wavelengths are detected as two components. The wavelengths of the lines are taken from the CXC atomic database 
ATOMDB and the table of Vemer et al. (1996), those of the series limits (= K-ionization thresholds) are from Vemer & Yakovlev (1995). 


TABLE 5 

Results from the detected absorption lines from H and He-like ions: velocity shifts v = (A-A 0 )/A 0 c in kms -1 


H-like 

O 

VIII 

Ne 

X 

Na 

XI 

Mg 

XII 

Al 

XIII 

Si 

XIV 

S 

XVI 

Ar 

XVIII 

Ca 

XX 

Fe 

XXVI 

Ly a 
Ly (3 
Ly 7 
Ly 5 
Lye 

71 q+1716 
/18 -281 

-90 +187 

~ -185 

—70+85 
J -81 

-109 +1109 

1UV -198 

-128+35 

79+50 

/Z -47 

7+64 

-64 

i Q9+46 
iyZ -46 

288+ 2 7 3 

284^1 

-34^ 

ti+922 

-1401 

— 171+88 
1 ' -104 

744+431 

J ^-448 

-601 2 

7C+263 
' -238 

308! 333 

-41+114 

-95 

-147+ 393 

^-399 

49Q+536 

^ zu -489 

-3221! 2 

-7501$ 

164 1o! 


Ly series 

-41 +6 ? 
Hi -72 

-97+28 

y/ -20 

207!}gt 

-38+H 

-183+}$ 

-54+33 

J -32 

-90+!! 

-327l‘ g 





o 

Ne 

Na 

Mg 

Al 

Si 

S 

Ar 

Ca 

Fe 

He-like 

VII 

IX 

X 

XI 

XII 

XIII 

XV 

XVII 

XIX 

XXV 

He a 


-203+ 3 ! 


_71+35 
' l -49 

286H 

- 301 ! 

84+ 127 
° -138 

_142+ 250 

-233 

1124 +98 
1 -1069 

_116+ 501 
1 1 -531 

He/3 

507 +492 

-1020 

-361 7 


-651 4 


1811 

609+l|| 6 

79+440 

-306 

—74+616 
' ^22 


He 7 


-507+!!! 


77+81 

' -71 


-19+320 
1 -234 

7021 7 ! 

655 +283 



He 5 











He e 


-Ql+54 
" -59 









He series 

-291 ±(>1000) 

-158^56 


-72+25 
' -28 

1351“ 

-86l« 


9+8 

-12 


-21 +0 

z,1 -379 


A negative velocity indicates a blue shift, due to the absorbing material moving towards the observer. Rows labeled with “Ly/He series” show the results from 
modeling the complete absorption line series of the corresponding ion at once with a single physical model (§3.5). 
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Fig. 10. — The Chandra (non-dip) spectrum of CygX-1, shown as ratio of data and continuum-model (Table 2). For visual clarity, the data have been rebinned 
to a common resolution of 10 mA, and all MEGdzl and HEGzbl spectra have been combined. The Gaussian line fits and identifications are shown as well. (The 
labels mark the lines’ rest wavelengths.) Iron L-shell transitions (lines of Fe <XXV) are shown in blue in the electronic version. 
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TABLE 6 

Results from the detected absorption lines from H and He-like ions: column densities in 10 l6 cm~ 2 


H-like 

O 

VIII 

Ne 

X 

Na 

XI 

Mg 

XII 

Al 

XIII 

Si 

XIV 

s 

XVI 

Ar 

XVIII 

Ca 

XX 

Fe 

XXVI 

Ly a 
Ly/3 

L y 7 

Ly S 

3 1 

151 

201 

H+J3 
1 -20 

3 5 30 ' 2 
tO.81; 2 

511 

82j“ 


5-21J 

8 °-i;f 

1.8 ±0.5 

3+3 

-2 

7 -i-81 
181 
19 ±13 

5 0+ 0 - 6 

J -l.i 
18-15 

29™ 

2 2 +1 - 2 
z>z - 1.8 

7-7 

2 ± 2 


Ly series 

36 +18 

-12 

42+ro 

9 q+0.4 
^• d -0.5 

7.2 ±0.6 

1.4 ±0.5 

10.1 ±0.8 

15 j* 2 

111 




He-like 

O 

VII 

Ne 

IX 

Na 

X 

Mg 

XI 

Al 

XII 

Si 

XIII 

s 

XV 

Ar 

XVII 

Ca 

XIX 

Fe 

XXV 

He a 
He/3 
He 7 

41 

1.5 ±0.2 

9 O+ 0.8 

2 '°-0.4 

3.0± 1.5 


1 57+ 0 - 13 
5.6 ±0.9 
8±2 

o-S81 

3 - 8 1:i 

i 04 +O.II 
1 -0.18 

9-3 

81 

9 9±0.4 

zo -0.5 

16jJ 

23 +n 
z -14 

1.3 ±0.6 
4±4 
21-1 

LOli 

111 

151 

He series 

01 

6 3 +1 - 3 
° -1.6 


5+2 

-1 


12 ±3 

15+ 2 



1461° 


The column densities for the single lines have been calculated using Eq. (11), assuming that the line is on the linear part of the curve of growth, which 
underestimates the column density for saturated lines. Rows labeled with “Ly/He series” show the results from modeling complete absorption line series (§3.5). 


be inferred (Spitzer 1978, Eq. 3-48; Mihalas 1978, §10-3): 


mc 2 -W x = 1.13x 10 17 cm~ 2 ( ( Wx\ 

fj \Aj\mAj 


7 T 


If the lines are, however, saturated, W\ depends on Aud as 
well, and one has to construct the full “curve of growth” with 
several lines from a common ground state i in order to con- 
strain Nj (see, e.g., Kotani et al. 2000). 

We have therefore performed a systematic analysis of ab- 
sorption line series: H-like ions are detected by their Is — > np 
Lyman series, He-like ions by their Is 2 — ■» Is up resonance ab- 
sorption series, see Table 4. For those lines which are clearly 
detected and not obviously affected by blends, the measured 
velocity shifts (Eq. 8) are shown in Table 5. Most of the lines 
are detected at rather low projected velocity ( | v| <200 km s -1 ). 
Note that the systemic velocity of Cyg X-l/HDE 226868 is 
(-7.0 ± 0.5) km s -1 (Gies et al. 2003), and that the radial ve- 
locity of both the supergiant and the black hole vanishes at 
orbital phase cf> OI b = 0, while the radial component of the fo- 
cused stream should be maximal (likely up to 720 km s -1 , see 
§4.1). The column densities in Table 6 are calculated using 
Eq. (11), assuming that the line is on the linear part of the 
curve of growth. As the strongest lines are, however, often 
saturated, Eq. (11) predicts too small column densities from 
them. For weak lines, on the other hand, the equivalent width 
is most likely overestimated such that the quoted values may 
rather be upper limits. The properties of the lines (Einstein-A 
coefficients and quantum multiplicities, which determine the 
oscillator strength, as well as the rest wavelengths) are taken 
from Verner et al. (1996) and ATOMDB 13 v. 1.3.1. 


3.5. Line series of H-/He-like and Fe L-shell ions 

As an alternative to the curve of growth, we chose to de- 
velop a model which implements the expected line profiles 
(Eq. 10) directly for all transitions of a series from a common 
ground state i. The model contains A/, £o, and the systemic 
shift velocity v (Eq. 8) as fit parameters, and thus avoids the 
use of equivalent widths at all. This approach allows for a 
systematic treatment of the iron L-shell transitions as well, 

13 See http://cxc.harvard.edu/atomdb/. 



Wavelength [A] 


Fig. 1 1 . — The absorption line series of Fe XIX between 12.8 A and 14 A. 
The Ne IX r line which blends with the complex at 13.4-13.6 A is also shown, 
and the expected positions of the Ne IX i and f emission lines are indicated. 


TABLE 7 

Parameters of the absorption line series for Fe L-shell ions 


ion v A(ion) £o 

(fans -1 ) (10 16 cnT 2 ) (kms~‘) 


Fe XXIV 

21 +31 

Zi -30 

A O+0.7 

82l> 

Fe XXIII 

12015 

9 9 ±0.4 
z ' z - 0.2 

41011 

FeXXII 

_27+ 33 

z -35 

9 q± 0.2 
z -0.4 

19311 

Fe XXI 

-156! 3 “ 

2 5+ 0 ’ 2 
z " -0.4 

26«1“ 

Fe XX 

-2831 3 

3 0 + 0 ’ 3 
j.u _0 2 

333l 2 

Fe XIX 

- 1161 “ 

3 1+ 0 - 2 
■^-oj 

251!§ 

Fe XVIII 

-60l> 

i o±0.3 
1 '°- 0.2 

20 3 1 | 

FeXVII 

-2771“ 

1.4±0.3 

78! 2 “ 


which are often blended with other lines such that the different 
contributions to a line’s equivalent width can hardly be sepa- 
rated when only single Gaussians are used. As an example. 
Fig. 1 1 shows the Fe XIX complex between 12.8 A and 14 A. 
Lines from different Fe XIX transitions overlap, and so does 
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9.15 9.2 9.25 9.3 9.35 


Wavelength [A] 

Fig. 12. — The Mg XI triplet of resonance (r), intercombination (i) and 
forbidden line (f) in the Chandra-HETGS spectrum. The MEG data is shown 
with an offset of —0.5. While the ion/line labels indicate the rest wavelengths, 
the abbreviations m, in and in label the actual line positions (see Table 8). 


TABLE 8 

Parameters of the lines of the Mg xi triplet 


line 

component 

A 

V 

Wx 



(A) 

(kms -1 ) 

(mA) 

Mg XI r 

l-i 

9 167+ 0001 
yAO 7 -0.002 



Mg XI r? 

(r2) 

Q 1QO+0.004 

(700!‘2°) 

-3.5 ±0.9 

Mg XI i 

il 

9.232 + 0.003 

U+107 
1 -105 

+1 5 +l5 
‘ A ■‘L-0.6 

Mg XI i 

i2 

9 261 +0 002 

5LZ,Ul_oo03 

968^o 

+1 7 +11 
+i,/ -0.9 

Mg XI f 

fl 

Q 1+0-004 

^ -0.001 

-110!* 30 

+2 5+ 0 - 9 
-t-z.j— 1 2 

Mg XI f 

f2 

Q QOC+0.009 
“•• z,z ' o _0 003 

454+2” 

+1 3^- 9 
-1.0 


The shift velocity v and equivalent width W\ are given by Eqs. (8) & (9). 
Note that the (r2) component is likely caused by absorption from blueshifted 
transitions of Fe XXI, Fe XX and Ne X. 


the strong Ne IX r-line. Furthermore, the absorption features 
are often rather weak and no prominent lines can be fitted, 
whereas the line series model can still be applied. Although 
a disadvantage of this approach is the larger computational 
effort, it usually allows to constrain the parameters of a line 
series more tightly. 

The model with physical absorption line series fits the data 
hardly worse than the model with single Gaussian lines: A \ 2 
of 12812 instead of 12180 before (see Table 2) is obtained. 
The results are presented in the last row for the whole Ly/He 
series of Tables 5 & 6 for the H-/He-like ions, and in Table 7 
for the Fe L-shell ions. The column densities inferred from 
the series model are generally in good agreement with the val- 
ues derived from the single Gaussian fits for the higher transi- 
tions of H- or He-like ions (Table 6), while the a - sometimes 
even (3 - lines are saturated. These measurements can be used 
as input for wind or photoionization models, although a few 
columns are rather badly constrained if the thermal velocity £o 
is left as a free parameter: As the line profiles are not resolved, 
there is a degeneracy between Doppler broadened lines and 
narrow, but saturated lines. The notable fact that no strong 
wavelength shifts are observed is, however, independent of 
this degeneracy: almost all line series are consistent with a 
velocity between -200 km s -1 and +200 km s _1 . Lower ion- 
ization stages of the same element are usually seen at higher 
blueshift, like in the sequence Fe XXIII - XXII - XXI - XX. 


TABLE 9 

Electron densities corresponding to R = F ( f )/ F ( i ) ratios 




model calculations by 
Porquet & Dubau (2000) fl 



R { Mg XI) 

n e 

component 

measurement 


(10 12 cm -3 ) 

first pair 


2.3 

1.0... 4.5 

of i & f lines 

R, (Mg XI) =1.6^ | 

1.6 

3.5... 10 

(unshifted) 


1.0 

8 ... 22 

second pair 

r 

1.3 

5... 15 

of i & f lines 

R 2 ( Mg XI) = 0.8 ±0.5 < 

0.8 

12... 28 

(redshifted) 


0.3 

43... 91 


a A temperature between 0.3 and 8 MK is assumed 
and UV -photoexcitation is not considered. 


3.6. Emission lines from He-like ions 

The transitions between the Is 2 ground state and the ls2s 
or ls2p excited states of Helium-like ions lead to the triplet of 
forbidden (f), intercombination (i), and resonance (r) line, see 
Table 4. These lines provide a density and temperature diag- 
nostics of an emitting plasma via the ratios R(n e ) := F( f) / F( i) 
and G(T e ) := [F(i)+/ r (f)] / F( r) of the fluxes F in the r-, i-, and 
f-line (see, e.g., Gabriel & Jordan 1969; Porquet & Dubau 
2000). In this observation, the dipole-allowed resonance tran- 
sitions are seen as absorption lines, as an absorbing plasma is 
detected in front of the X-ray source. For the same reason, we 
cannot use the Fe L-shell density diagnostics (Mauche et al. 
2005), as, e.g., the Fe XXII emission lines at 11.77 A and 
1 1 .92 A used by Mauche et al. (2003) are both seen in absorp- 
tion. But we can still use the detected He-i and -f emission 
lines to estimate the density via the /Gratio, noting the caveat 
that the densities are systematically overestimated in the pres- 
ence of an external UV radiation field, as photoexcitation of 
the ls2s ( 3 Si) — > ls2p ( 3 P) transition depopulates the upper 
level of the f-line in favor of the i-line, and leads to a lower 
/Gratio (see, e.g., Mewe & Schrijver 1978; Kahn et al. 2001; 
Wojdowski et al. 2008). 

The i- and f-line of the Mg XI triplet are seen as two dis- 
tinct components each - one with almost no shift, which is 
consistent with the r-absorption line, and the other one at a 
redshift of 400- 1000 km s _I , see Fig. 12 and Table 8. Given 
these two emission components and their Doppler shifts, 
one could be tempted to identify the absorption feature at 
9. 19 A with a second redshifted MgXlr line, but our model 
for the complete absorption line series (§3.5) predicts that 
the blueshifted 2s — > 4p transitions from the ground states 
of FeXXl and FeXX, as well as the Nex 1 s — > 1 Op transi- 
tion, account for most of the absorption seen at 9.19 A. Ta- 
ble 9 shows the R -ratios obtained for the two pairs of lines 
as well as the corresponding densities according to the cal- 
culations of Porquet & Dubau (2000, Fig. 8e), neglecting the 
influence of the UV radiation of HDE 226868 on the R-ratio. 
The unshifted lines would then be caused by a plasma with 
an electron density n e , l ~ 5 x 10 12 cm -3 , the redshifted pair 
of lines seems to stem from another plasma component with 
« e . 2 « e ,i- A more detailed discussion of these densities is 

presented in §4.1. 

The S/N ratio of the spectrum above 20 A is not high 
enough to describe the O VII triplet. The Ne IX triplet blends 
with several Fexix lines (Fig. 11). Nax i- and f-lines 
are likely present with several components, but those can 
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not clearly be distinguished. Similar as for Mg, there are 
also two i-lines of Alxil at shifts of -477^40 km s _1 and 
+333^*28 kms -1 , respectively, but no f-line is detected due to 
a blend with the strong Mg XI He (3 absorption line and es- 
pecially a Fe XXII ls 2 2s2p(2s^5p) absorption line at 7.87 A. 
Similarly, the Si XIII i-line is not detected as it overlaps with 
absorption features probably due to nearly neutral Si K-edge 
structures. (Furthermore, the flux of the Si XIII f-line might be 
underestimated as it blends with the Mg XII Ly7 line.) The 
i- and f-lines of He-like sulfur are detected with R ( S XV) = 
0.61^ ^, but S XV has not been modeled by Porquet & Dubau 
(2000). Below 4 A, the S/N ratio below is again not high 
enough to resolve further triplets from heavier elements such 
as argon, calcium or iron. 

4. DISCUSSION & CONCLUSIONS 

In the following, we discuss our results and derive con- 
straints on the stellar wind in the accretion region. 

4.1. Velocity and density of the stellar wind 

While a spherically symmetric model for the stellar wind 
in the HDE226868/CygX-l system can be excluded by ob- 
servations (see, e.g., Gies & Bolton 1986b; Gies et al. 2003; 
Miller et al. 2005; Gies et al. 2008), a symmetric velocity law 

v(r) = Voo ■ f(r/Rf) (12) 

is usually assumed to obtain a first estimate of the particle 
density in the wind. The fraction / of the terminal velocity v'oo 
is often parameterized by (Lamers & Leitherer 1993, Eq. 3) 

f(x) = f 0 + (l-/o)-(l-l/x) /3 (for x > 1) (13) 

with /o := vo/vqo, where vo is the velocity at the base of the 
wind (x = 1). The simple model for the radiatively driven wind 
of isolated stars by Castor et al. (1975) is obtained for (3= 1/2. 
Photoionization of the wind, however, suppresses its acceler- 
ation (Blondin 1994), such that a smaller / (e.g., a larger (3 
within the same model) is required in the Stromgren zone. 
Gies & Bolton (1986b) have explained the orbital variation of 
the 4686 A Hell emission line profile of CygX-1 with a sim- 
ilar model for the focused wind, where Voo, /?*, and f3 depend 
on the angle 6 from the binary axis. (3 was interpolated be- 
tween 1.60 and 1.05 for 8 between 0 and 20°. The value (3(6= 
20° ) = 1.05 is, however, often used for a spherically symmet- 
ric wind as well (Lachowicz et al. 2006; Szostek & Zdziarski 
2007; Poutanen et al. 2008). Vrtileketal. (2008) use /o = 
0.01 to avoid numerical singularities and fit the (relatively 
low) value f3 ~ 0.75 to their models for UV lines, vo is 
likely of the order of the thermal velocity of H-atoms, which 
is {2kT = 23 kms -1 and corresponds to f 0 = 0.011 
for HDE226868’s effective temperature T e g- = 32000K and 
Voo = 2100 kms -1 (Herrero et al. 1995). 

With the mean molecular weight per H atom, /.t « 1.4, 14 
the continuity equation M* = / j,mnnn(r ) • 4-7rr 2 v(r) gives the 
following estimate for the hydrogen density profile: 

"h (r) = -d{r/Rf) with d{x) = x~ 2 /f(x) (14) 

4ttR^Voo 

Using the parameters of HDE 226868 (Af+ = 3 x 10 -6 Mq yr -1 , 
= 11 Rq\ Herrero et al. 1995; also summarized by 
Nowak et al. 1999, Table 1), Eq. (14) predicts 

n H (r) = 2.2xl0 10 cm -3 ■ dfr/Rf) . (15) 

14 The helium abundance per H-atom is m\0% (Wilms et al. 2000). 


TABLE 10 

Solutions to Eq. (15) for the model ofEq. (13) 


«h(x) 

(10 10 cm -3 ) 

d(x) 

A 

= Vq/Voo 

p 

II 

* 

>3 

2100 -fix) 
for j3 = /3max 

440 

200 

0.005 

< oo 

l 

11 

110 

50 

0.011 

< 1 

< 1.01 

<41 

no 

50 

0.011 

<2 

< 1.08 

<36 

110 

50 

0.011 

< 3 

< 1.18 

< 30 

22 

10 

0.011 

< 1 

< 1.08 

< 180 

22 

10 

0.011 

< 2 

< 1.29 

< 127 

22 

10 

0.011 

< 3 

< 1.48 

< 95 

4.4 

2 

0.011 

< 1 

< 1.36 

< 570 

4.4 

2 

0.011 

< 2 

< 1.69 

< 368 

4.4 

2 

0.011 

< 3 

< 1.97 

<271 


Note that the binary separation a = 41 Rq corresponds to x a = a/R* = 2.4. 


Eq. (14) can only be solved within the model of Eq. (13) 
- or any other model for / in Eq. (12) with /(x) > /o - for 

d < Jo 1 and 1 < x < y/l/(fod) . (16) 

For /o«0.01, the value d « 1 90 - which would be required to 
explain the density 15 «h.i ~ 4.2x 10 12 cm -3 obtained from the 
unshifted Mg XI triplet (§3.6) - can never be reached within 
our model for the continuous spherical wind. This result 
shows that the density is likely overestimated by an /C ratio 
analysis which ignores the strong UV-flux of the 09.7 star. 
Table 10 lists therefore some solutions to Eq. (15) for much 
lower densities as well. Due to the additional constraint that 
the radial velocity is < 100 km s -1 (Table 8), we suggest that 
the first emission component stems from close to the stellar 
surface. Although the simple wind model of Eq. (13) may not 
be appropriate in this region, the results of Table 10 are rather 
insensitive to the assumptions on the velocity law, as a wide 
range of (3 values was considered, but depend mostly on the 
wind’s initial veloctiy vo, for which we have used a reason- 
able estimate. Note that /o and d depend also on Voo, but their 
product, which is important in Eq. (16), does not. 

In spite of the systematic errors of the (absolute) density 
analysis with the /Uratio, we infer that the second plasma 
component is much denser relative to the first one. As it 
is seen at a larger redshift of 400- 1000 kms -1 , we favor 
its identification with the focused wind. The two emission 
components could also be caused coincidentally by dense 
clumps in the stellar wind (which are common for O-stars; 
e.g., Oskinova et al. 2007; Lepine & Moffat 2008), but the 
interpretation as a slow base of the wind close to the stel- 
lar surface and a focused wind between the accreting black 
hole and its donor star provides a consistent descripion of 
both emission components: An undisturbed wind (with v x 
and /o as above, and (3 = 1 .05 — > 1 .6) would reach a veloc- 
ity of via) « 900 <— 1200 km s -1 in the distance of the black 
hole. The focused wind in the orbital plane would then be de- 
tected with a projected velocity v(a)sin; = 500 <— 700 kms -1 
at </> or b = 0. While photoionization reduces the effciency of 
acceleration for the spherical wind, the denser focused wind 
is less strongly affected due to self-shielding and reaches the 
expected velocity. It is an observational fact that the fo- 
cused wind has another ionization structure: Optical emis- 
sion lines (Ha, He II A4686) from ions which only exist at a 
low ionization parameter have been observed in the focused 

15 A plasma of fully ionized H and He contains -e 1 .2 electrons per H-atom. 
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Fig. 13. — Relative population of the ionization stages of iron as a function 
of the ionization parameter £ = L/(nyir 2 ), as calculated with XSTAR for 
L 37 = 3.5 and = 0.01 . 

stream (see, e.g., Gies & Bolton 1986b; Ninkov et al. 1987; 
Gies et al. 2003). For the spherical wind, however, Gies et al. 
(2008) have conjectured that the part between CygX-1 and 
the donor star might be completely ionized in the soft state. 
Our observations show that the situation in the hard state may 
be similar. 

4.2. Modeling of the photoionization zone 

We use the photoionization code XSTAR 2.11n7b 16 
(Kallman & Bautista 2001) to model the photoionization 
zone. As the latter is quite complex (due to the inhomoge- 
neous wind density, which is strongly entangled with the X- 
ray flux), we do not claim to describe the photoionized wind 
self-consistently, but only want to derive a first approxima- 
tion. For an optically thin plasma, the relative population of 
a given atom’s ions is merely a function of the ionization pa- 
rameter 

= Li \ erg cm s -1 , (17) 

n H r z «13G 2 

where L37 is the ionizing source luminosity above 13.6eV 
in 1 0 37 erg s -1 , « 13 is the hydrogen density in 10 I3 cm -3 and 
ri2 is the distance from the source in 10 12 cm. 

For Cyg X-l as observed in the hard state, extrapolation of 
the unabsorbed model obtained by our analysis (Table 2) gives 
L37 ~ 3.5. Although there are obviously strong variations 
in the density, the XSTAR calculation has been performed 
with constant = 0.01, which is the average of /; H (rj from 
Eq. (15) for R+ < r < a. Fig. 13 shows the resulting popula- 
tion of Fe ions. The ionization parameter at which the pop- 
ulation distributions of some ions of interest peak, as well as 
the corresponding full widths at half maximum (FWHM) are 
presented in Table 11. We have also included Si IV, N V and 
C IV for a comparison with the work of Vrtilek et al. (2008) 
and Gies et al. (2008). The H- and He-like ions detected with 
Chandra only appear at considerable distance from the X-ray 
source. Taking into account that the photons emanating from 
Cyg X-l first propagate through a lower density wind and that 
the (eventually stalled) wind of higher density is only reached 
in the vicinity of the star, the actual distances will even be 
larger than the r\i values quoted in Table 1 1 . 

16 See http://heasarc. gsfc.nasa.gov/lheasoft/xstar/xstar.html. 


TABLE 1 1 

Ionization parameter for peak ion populations 


Ion 

log? 

FWHMflogO 

ri2=(G7/«13/10‘°8«) 1/2 

Fe XXVI 

2.91 

0.27 

0 . 6 - 0.8 

Fe XXV 

2.75 

0.32 

0.7-1 .0 

FeXXIV 

2.65 

0.30 

0 . 8 - 1. 1 

Fe XXIII 

2.58 

0.28 

0 . 8 - 1. 1 

Fe XXII 

2.51 

0.31 

0.9-1. 3 

Fe XXI 

2.44 

0.36 

0.9-1 .4 

Fe XX 

2.32 

0.44 

1.0-1 .7 

Fe XIX 

2.12 

0.40 

1.2-1 .9 

Fe XVIII 

2.01 

0.23 

1 . 6 - 2.0 

Fe XVII 

1.97 

0.14 

1 . 8 - 2. 1 

CaXX 

2.69 

0.56 

0.7-1 .2 

CaXIX 

2.39 

0.54 

0.9-1 .6 

ArXVIll 

2.55 

0.65 

0.7-1 .5 

ArXVII 

2.22 

0.50 

1 . 0 - 1. 8 

Sxvi 

2.38 

0.69 

0.8-1 .7 

S XV 

2.08 

0.43 

1 . 2 - 2.0 

Si XIV 

2.21 

0.65 

0.9-2.0 

Si XI 1 1 

1.96 

0.39 

1. 5-2.4 

Si IV 

1.09 

0.18 

4.7-5.7 

Mg XII 

2.06 

0.58 

1. 2-2.3 

Mg XI 

1.82 

0.31 

1 . 8 - 2.6 

NeX 

1.94 

0.40 

1. 6-2.5 

NeIX 

1.74 

0.24 

2 . 1 - 2.8 

O VIII 

1.80 

0.26 

1 . 9-2.6 

O VII 

1.62 

0.22 

2.5-3.2 

Nv 

1.51 

0.01 

3.25-3.30 

Civ 

1.48 

0.19 

3.3-4. 1 


We list the range in ionization parameter and distance for the maximum 
ion population (Fig 13) obtained in an XSTAR simulation with L 37 = 3.5 
and n 13 = 0.01. Note that the binary separation in units of 10 12 cm is ai 2 = 2.9. 

4.3. Origin of redshifted X-ray absorption lines 

Additional constraints on the accretion flow can be derived 
by considering the Doppler shifts observed in absorption lines 
(§3.4 & 3.5). Models for the wind velocity like the one of 
Eqs. (12) & (13) predict a velocity v(r) at the position r in the 
stellar wind. But only the projection of v against the black 
hole, v ra d,BH, can eventually be observed as radial velocity in 
an X-ray absorption line. With the angle a{r ) between v(r) 
and the direction from r towards the black hole, v ra d.Bii(r) is 

i'rad,BH(r) = cos a{r) ■ |v(r)| . (18) 

Assuming a velocity field with radial symmetry with respect 
to the star, a(r ) is just the angle at r between the star and the 
black hole. For example, a is 90° (and v r ad,BH is therefore 0) 
on the sphere containing the center of the star and the black 
hole diametrically opposed. Redshifted absorption lines can 
only be observed from wind material inside this sphere, while 
the part of the wind outside of it is always seen at a blueshift. 
- a fact which is independent of the assumed velocity field as 
long as it is directed radially away from the star. 

For a spherically symmetric wind model, any line of sight 
can be rotated to an equivalent one in a half-plane limited 
by the binary axis. Only a sector with half-opening angle i 
can be observed unless the system’s inclination is i = 90°. 
Fig. 14 shows the projected velocity vVad.Bibf) for a simple 
wind model and two sectors containing all observable lines of 
sight towards CygX-1 for an inclination of i = 35°. 

We now investigate the region where the projected wind ve- 
locity (Fig. 14) is compatible with the observed Doppler shifts 
(Table 5 and Table 7). We are confident that this method al- 
lows for the identification of the absorption regions, as the 


14 


M. Hanke et al. 



Fig. 14. — Gray scale image of v r ad BH(r) of Eq. (18), i.e., projected 
wind velocity against the black hole. A spherically symmetric velocity law 
given by Eqs. (12) & (13) with Voo = 2 1 00 km s f fy = 0.01 and (3 = 1.05 is 
assumed. The star and the black hole are shown as filled black circles; the size 
of the latter is the accretion radius 2GMbh / v(a) 2 . On the black circle passing 
through the center of the star and through the black hole v rat ] bhM is 0 since 
a(r) = 90°. Positive velocities (+) are seen only within this circle', a lighter 
gray means a larger redshift. Negative velocities (— ) can likewise only occur 
outside of this circle', a lighter gray here means a larger blueshift. Dashed gray 
lines of constant v ra d,BH are shown from -1400km s 1 to +1000 km s 1 in 
steps of 200 km s -1 . The two highlighted sectors contain all observable lines 
of sight towards the black hole (after rotation in this plane) for the inclination 
i = 35° . The labels show the corresponding orbital phases. 

low observed velocities are always found close to the a = 90° 
sphere - independent of the wind model. For most of the in- 
vestigated ions (Ne x. Mg xn. Si xiv & xin, S xv, Ar xvii, 
Fe XIX & XVIII) the inferred distance from the black hole 
agrees with the predictions of Table 11; e.g., the projected 
velocity -117 kms -1 < VNeX <-69 kms -1 measured for Ne X 
is (during 0.93 < </> or b < 1) obtained at m = 1.78 ±0.07. 
For many other ions, both results are still very similar; e.g., 
Neix with -214kms _1 < VNeix < -123 kms -1 is expected at 
r\2 = 1.95 ± 0.07 from the v ra d,BH-model and at r\s = 2.1 -2.8 
from the photoionization model. Only for the highly ionized 
iron lines, the small distances are - as already anticipated in 
§4.2 - underestimated by the XSTAR simulation run with 
constant average density, which overestimates the wind den- 
sity close to the black hole. For example, the velocity range 
-9km s -1 < VFexxiv < 52km s -1 measured for Fe XXIV corre- 
sponds to r \2 =1.51 ±0.07, while the population of this ion 
peaks at = 0. 8-1.1 within the model presented in Table 1 1 . 


4.4. Previous Chandra observations 

In previous Chandra -FIETGS observations of Cyg X-l, the 
stellar wind was seen under different viewing angles, or the 
X-ray source was in other spectral states (see Fig. 15), which 
probably means that the properties of the wind were also dif- 
ferent. Schulz et al. (2002) have analyzed the 14 ks obser- 
vation #107 performed in 1999 October (at (j> or b « 0.74), 17 
when Cyg X-l was in a transitional state. They derive a neu- 
tral column density of Ah = 6.2 x 10 21 cm -2 from prominent 

17 Note that Schulz et al. (2002) and Miller et al. (2005) quote an erro- 
neous date and thus the wrong orbital phase <f > or b - 0.93 f 0.74 for this obser- 
vation. 
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Fig. 15. — Orbital phase and mean ASM count rate of the previously re- 
ported Chandra observations of CygX-1. (#3814 is presented in this paper). 

absorption edges (see also Table 3) and detect some emis- 
sion and absorption lines with indications of P Cygni pro- 
files. Miller et al. (2005) have investigated the focused wind 
with the 32 ks observation #2415 of CygX-1 in an interme- 
diate state, performed in 2001 January (at (/> OI t ~ 0.77). They 
report absorption and emission lines of H- and He-like res- 
onance lines of Ne, Na, Mg, and Si with a mean redshift 
of ss 100 kms -1 , as well as some lines of highly ionized Fe 
and Ni. The column density is 6.2 x 10 21 cm -2 (Miller et al. 
2002). Marshall et al. (2001) describe Ly a and He a absorp- 
tion lines, redshifted by (450 ± 1 50) km s~ from the 12.6 ks 
observation #1511 of CygX-1 in the hard state, performed 
in 2001 January (at (/) or b ~ 0.83). Chang & Cui (2007) report 
dramatic variability in the 30 ks observation # 3407 performed 
in 2001 October (at <j> or b ~ 0.88), when CygX-1 reached its 
soft state: While a large number of absorption lines (mostly 
redshifted, but not by a consistent velocity) is identified in the 
first part of their observation, most of them weaken signif- 
icantly or cannot be detected at all during the second part. 
Complete ionization of the wind due to a sudden density 
decrement is given as a possible explanation. Feng et al. 
(2003) detect asymmetric absorption lines with the 26 ks ob- 
servation #3724 of CygX-1 in the soft state, performed in 
2002 July (at <j> or b = 0). The line centers are almost at their rest 
wavelengths, but the red wings are more extended, especially 
for the transitions of highest ionized ions, which they explain 
by the inflowing focused wind reaching both the highest red- 
shift and ionization parameter closest to the black hole. 

The interpretation of the absorption lines presented in the 
previous section describes our observation consistent with 
wind and photoionization models. It can, however, not be 
applied to all of the other observations; the model of Fig. 14 
neither predicts a conspicuously high redshift at (f > or b ~ 0.83 
and at a higher soft X-ray flux (observation #1511), nor a pos- 
itive redshift at <p OI \, « 0.77 at a still higher flux (observation 
#2415, see Fig. 15). Inhomogeneities (e.g., density enhance- 
ments in shielding clumps) and asymmetries of the wind (e.g., 
due to non-inertial forces in the binary system) may therefore 
play an important role in these cases. 

4.5. Summary 

In this paper, we have presented a Chandra-U F I GS obser- 
vation of CygX-1 during superior conjunction of the black 
hole, which allows us to detect the X-ray absorption signa- 
tures of the stellar wind of HDE 226868. The light curve near 
Gorb = 0 is shaped by absorption dips; these are, however, ex- 
cluded from this analysis by selecting non-dip times of high 
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count rate only. At a flux of ~0.25 Crab, we have to deal with 
moderate pile-up in the grating spectra, for which can be ac- 
counted very well with the simple_gpile2 model. The 
continuum of both Chandra ' s soft and RXTE's broad band 
X-ray spectrum has been described consistently by a single 
model, consisting of an empirical broken power law spectrum 
with high energy cutoff (which is typical for the hard state) 
and a subordinated disk component with an inner radius close 
to the ISCO. The joint modeling of both spectra reveals the 
presence of both a narrow and a relativistically broadened flu- 
orescence Fe Ka line. Chandra has resolved absorption edges 
of neutral atoms with an overabundance of metals, suggesting 
an origin not only in the ISM, but also in the stellar wind of the 
evolved supergiant. The previously suspected anti-correlation 
of neutral column densities and X-ray flux, which is due to 
photoionization, is confirmed. Absorption lines of highly ion- 
ized ions are produced where the wind becomes extremely 
photoionized. For H- and He-like ions, Lyman and He-series 
are detected up to the Ly or He e line, which we use to measure 
the column density and velocity of absorbing ions. For the 
wealth of Fe L-shell transitions, column densities can best be 
obtained by directly using a physical model for the complete 
line series of an ion. The non-dip spectrum shows almost no 
Doppler shifts, probably indicating that we have excluded the 
focused wind by our selection of non-dip times. We have also 
detected two plasma components in emission by two pairs of 
i- and f-lines of He-like Mg XI. The first one is roughly at 
rest and we identify it with the base of the spherical wind 
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close to the stellar surface. The second plasma component 
is denser and observed at a redshift that is compatible with 
the focused wind. A simple XSTAR simulation indicates 
that most of the observed ions only exist in a distance of the 
black hole, where the velocity of a spherical wind, projected 
against the black hole, is low - which is a consistent explana- 
tion of the small Doppler shifts observed in absorption lines. 
We review the previously reported Chandra observations of 
Cyg X-l and find that not all of them can be described in the 
same picture, as the wind may be affected by asymetries or 
inhomogeneities. A detailed spectroscopy of the absorption 
dips, which might shed more light on the focused wind, will 
be presented in a subsequent paper. 
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